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We examine the transformation of noise under a quantum error correcting code (QECC) con- 
catenated repeatedly with itself, by analyzing the effects of a quantum channel after each level of 
concatenation using recovery operators that are optimally adapted to use error syndrome infor- 
mation from the previous levels of the code. We use the Shannon entropy of these channels to 
estimate the thresholds of correctable noise for QECCs and find considerable improvements under 
this adaptive concatenation. Similar methods could be used to increase quantum fault tolerant 
thresholds. 



Decoherence is a significant problem for implementa- 
tion of quantum computing. Quantum error correction 
P, 0] mitigates decoherence and gate errors by intro- 
ducing redundancy. In particular, for many codes, one 
qubit is encoded into n qubits, and each of those encoded 
into n qubits, so that after j levels of concatenation, n J 
qubits are used. In previous work, concatenated quantum 
codes were analyzed using a channel map formalism to 
evaluate their performance and to determine the thresh- 
old error probability for perfect fidelity in the infinite- 
concatenation limit 3, 4]. In order to optimize the recov- 
ery operators at each level of the code, we develop here 
an adaptive concatenation approach, which yields sig- 
nificantly larger values for the thresholds of correctable 
noise. This work extends previous work [5J, by formulat- 
ing it in the channel superoperator representation, which 
allows for arbitrary quantum channels. We introduce the 
idea of using the Shannon entropy of ensembles of chan- 
nels to calculate the thresholds of correctable noise; we 
use this to find thresholds for two quantum error correct- 
ing codes under two different types of noise. 



channels are trace-preserving and hermitian-preserving, 
this implies that one-qubit channels written in the Pauli 
basis have real entries and the first row is A/} CT = 8i a . 
On one-qubit, this is 
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A quasi-density matrix p is an unnormalized density 
matrix. It can be normalized to a density matrix ^-p, 
where p s = tr p = cj. Just as a density matrix can be 
written as a sum of quasi-density matrices, a channel can 
be written as a sum of quasi-channels [4| . Quasi-channcls 
map density matrices to quasi-density matrices, and 
therefore do not have to preserve trace. If A/" is a quasi- 
channel and p is a density matrix, then the probability of 
the quasi-channel occurring is p c = tr(A/p) = J2a NiaCc, 
which depends on the state p unless J\fj a = p c Sj y(T1 in 
which case it gives the channel -^-Af — j^j AT with prob- 
ability p c . 



I. BASIC NOTATION 

The n qubit Pauli operators are V n = {I, X, Y, Z}® n , 
We use the notation that XY = X <g> Y. Each pair of 
Pauli operators a, a' either commutes or anticommutes, 
where a G V n . We let the quantity T](a, a') be 1 if they 
commute, —1 if they anticommute. 

We write a density matrix on n qubits as p = 
X^o-p-p c a a - A channel is a map from density ma- 
trices to density matrices, and can be written as %{p) = 
J2i AipA\, where the Ai are Kraus operators that act on 
the state. This channel has the matrix representation 
N~ = J2i O(Ai), where 0{M) = M ® M^ T and AT is a 
matrix that acts on the vector form of p in the standard 
basis. By looking at how the 4™ Pauli operators a G V n 
perform under the map $(er), the channel can be written 
as a 4™ by 4™ superoperator in the Pauli basis. Because 
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II. OVERVIEW OF STABILIZER CODES 

An [[n, k,d]] quantum stabilizer code encodes k logi- 
cal qubits into n physical qubits, and has a distance d. 
Typically, a quantum stabilizer code is given in terms of 
n — k generators gi, and 4 fe encoded Pauli operators a. 

The gi generate the stabilizer group S, which is iso- 
morphic to Z^~ k . The k logical qubits encode into the 
2 fe -dimensional code space Cs, which has the property 
that if |r/>) G C s and s G S, then s\ip) = |^>. The 
elements of the set C(S) C V n are the 2 n+k Pauli opera- 
tors which commute with S. They send states in C's to 
states in Cs- Each of the 4 fc different equivalence classes 
of C(S)/S corresponds to one of the 4 fc encoded logical 
Pauli operators a that act on the logical qubits in the 
code space. 

To perform error detection, we measure each of the 
n — k generators <?;, projecting into either the +1 or 
— 1 eigenspace for each generator g^. This gives us a 
syndrome bit fli, which is if the state is in the +1 
eigenspace, and 1 if it is in the —1 eigenspace. These 
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Pi form the 2™ _fc possible error syndromes (3. Pq = = 
(0, 0, ... , 0) is the syndrome corresponding to no error. 

To recover from an error characterized by syndrome 
P, a recovery operator R(f3) must be chosen that returns 
the qubit to the code space. Since giR(P)\4>) = R(P)\if>), 
then rj(R((3),g l ) = Pi for all i. The 2 n+k possible choices 
for R(P) are all in the same equivalence class of 'P n /C(S). 



III. CHANNEL MAP FORMALISM 

We define a 4™ by 4 fe encoding operator £ that maps 
from a density matrix on k qubits to one on n qubits as 



1 



i—k 



:£p. 



This projects us into the code space of the code. In the 
Pauli basis, the 4 fc columns of 2 „ 1 _ fc £ act as the logical 
Pauli operators in the code space, and outside the code 
space. They arc defined as 



£i = + 9i) and £t = c£j. 



(1) 



Now that the logical qubits are encoded into the code 
space, assume that some noise TV acts on the full 2™- 
dimensional space of the n physical qubits of the code. 
The c/i operators are measured, giving an error syndrome 
p. Some unitary recovery operator R(P) is then chosen. 
We can write the error-correction process as 

where 0{R{P)) is the error recovery superoperator, Vp 
is a projection superoperator into the syndrome p space, 
and Vq = 2 „- k £ o £ T (transpose of £ followed by £) is a 
projection into the code space. £* decodes back into the 
original k qubit space. The whole process can then be 
written as 
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-£* oJ20(R(P)) oMo£. 



(2) 



This represents the noise on the logical qubits Q in terms 
of the noise on the physical qubits N '. We can decompose 
Q as a sum of contributions from each error syndrome: 







G u = 
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: £* oO{U) oMo£. (3) 



from the physical one-qubit noise on each qubit to the 
logical one-qubit noise on the encoded code space. This 
is useful for analyzing a code concatenated with itself an 
arbitrary number of times. A logical qubit on one level 
of the code is treated as a physical qubit at the next of 
the code, allowing us to concatenate a code with itself 
many times. In particular, the code will correct noise 
if ]im k ^ oo n c ° k (N'W) = / 4]. However, this assumes 
that the recovery operators do not depend on syndrome 
information from the previous levels of error correction. 
Furthermore, it does not give as high thresholds as the 
optimal adaptive concatenation described below. 



V. DIAGONAL NOISE 

A Pauli matrix a has the diagonal superoperator 0(a), 
where 0(a) a > a > = rj(a,a'). Suppose the noise is diagonal 
in the Pauli basis. Then the Q u from Eq. [3] become 
Q%/ a , = 2^fc££' ° C( cr ) ° N ° Since the nonzero 

entries of the column £ CT ' are a's where s £ S, it follows 
that 



on— k / 4 
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a' 8,<j' s 



(4) 



If a quasi-channcl is diagonal, it can be written as a 
sum of Pauli superoperators JV = VaO{a\ which 
corresponds to having each Pauli error a with proba- 
bility p a . For a one-qubit quasi-channel with probability 
p = ^2„Pcr, the diagonal parts are 



\p,x,y,z] =pi[l, 1,1,1] +px[l, 1,-1,-1] 

+ P y[l,-l,l,-l],+p Z [l, -1,-1,1]. 

This yields the Pauli probabilities: 



Pi 
Py 



p + x + y + z 



p — x + y — z 



Px = 



(5) 



(6) 



Pz 



p — x — y + z 



Since the information entropy (Shannon) in terms 
of Pauli errors with probabilities qi is h{q{) where 
h(x) — —x\og 2 x, a quasi-channel gives a contribution 
of PY.a H'-y) = ~Hp) + J2a KPt) t0 the entropy of a 
channel. 



VI. ADAPTIVE CONCATENATION 



IV. SAME NOISE ON EACH QUBIT 

Often it is assumed that each qubit experiences the 
same noise, so N = (A/^ 1 -') 8 " '. Previous work 0, 0] 
showed that the case of a code C where k = 1 produces 
a map OP : A/" (1) -> £/3 ^r£* o 0{R{p)) o (,yV^)® n o£ 



We now show how using syndrome information from 
the lowest levels of a code can be used to optimize the 
threshold of the code. 

When we use the syndrome information, instead of 
having the logical encoded channel map written as a 
simple sum of channel contributions Q = Q R {(3) as 
in Eq. [31 we represent this additional information as 
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Y.0 G R{0) ® P- If we perform error correction and ig- 
nore the syndrome information, we have simply Q. The 
probability of measuring a syndrome (3 depends on the 
density matrix state p, and is 

^=tr(^)p) = 2«^^ ) c <7 . 

If the first row of this quasi-channcl in the Pauli basis 
is all zero except for the first term, i.e., Gf£ — pSi a , 
then pp = Qff = p. If the syndrome (3 is measured, the 
resulting noise is —Q R ^\ 

VII. RECOVERY OPTIMIZATION 

Given an [[n, 1, d]] quantum code, assume that the syn- 
drome (3 is measured. It has a corresponding equivalence 
class V n /C{S). From this, we choose some representative 
element rp to return to the code space. This element dif- 
fers from the optimal recovery operator R(f3) by a logical 
Pauli operator a. The recovery operators for (3 are there- 
fore equivalent to the 4 recovery operators &rp, where a 
is any one of the 4 logical encoded Pauli operators. We 
then have the quasi-channel Q Tfl . 

Now, for each of the n blocks of the code, we pass this 
quasi-channel Mi onto the next level of the code, and 
have the resulting noise TV = §Q;Mi. To optimize the 
recovery operator based upon this noise, either we can 
find the optimal arp at each level of the code, or we can 
apply a random recovery operator rp at each level of the 
code, and optimize at the end. In the second case, since 
each block differed by a Pauli operator <ii from optimal, 
the resulting noise differs by ®j0j from optimal. This 
is equivalent to having the recovery operator ®j<7j o rp 
instead of rp. This is aor a for some logical Pauli operator 
a and some error syndrome a, and therefore differs by a 
logical Pauli operator from the optimal recovery operator 
at the top level of the code. 

While the first method may seem more intuitive, the 
second may be more computationally efficient, since de- 
termining the optimal recovery operator at any given 
level of the code is nontrivial. At the last level of the 
code, we will choose the encoded Pauli operator a that 
yields the channel closest to the identity channel, which 
corresponds to the maximal p a from Eq. [6] If the re- 
sulting noise is very close to the identity channel, then 
we have established that very little logical noise is intro- 
duced at each step of a quantum computation, allowing 
long quantum computations to be performed. In prac- 
tice, the noise at the end of the concatenation will often 
be close to diagonal even though the initial physical qubit 
noise was not. 

Suppose we wish to calculate the exact optimized chan- 
nel map of an [[n, 1, d]] quantum code after j levels of con- 
catenation with itself. There are 2™ J_1 syndromes. Many 
of these syndromes may have the same channel compo- 
nents Q R<< P\ rendering it unnecessary to calculate all of 



TABLE I: Encoded Pauli operators a and encoding operator 
E a for the 2 qubit code. 



a 


a Scr 


I 

X 
Y 
Z 


II II + ZZ 
XX XX - YY 
XY XY + YX 

IZ IZ + ZI 


TABLE II: Q H W for diagonal noise, using the shorthand a = 
N a ,„. 


f3 R{(3) 




II 

1 XI 
1 IX 


| [II + ZZ, XX + YY, XY + YX, IZ + ZI] 
| [II - ZZ, XX - YY, XY - YX, IZ - ZI] 
i [II - ZZ, XX - YY, YX - XY, ZI - IZ] 


them, reducing 


the complexity. However, for many codes 



it is still computationally too expensive to compute the 
optimized channel past the second level. Nevertheless, 
the overhead involved with implementing the optimiza- 
tion is low because the encoded channel from one level 
is simply passed on to the qubits of the next level of the 
code. This leads to efficient Monte Carlo simulation. 



VIII. EXAMPLE: THE TWO QUBIT BIT FLIP 
CODE 

It is useful to look first at the two qubit bit flip code 
because of its simplicity; it is the classical two bit rep- 
etition code. Its code space is spanned by the logical 
encoded states |0) = 1 00) , and |1) = |11). The stabilizer 
group S = {II, ZZ} is generated by the one generator 
cji = ZZ. Tab. [T] lists the encoded Pauli operators and 
corresponding £ a (Eq. [1} . Using Eq. [31 we can calculate 
various channel map components G R ^K For example, 
Sx X z = \{-Nxx,iz + N XXtZI + N YYJZ - N YY . ZI ). In 
the case of diagonal noise, the channel map components 
can be found directly from Eq. 0J These are given in 
Tab. HH using the shorthand M a , a = a - 

Suppose now that we have the same diagonal noise 
M^ — [1, 1, X, x] on each qubit, with x > 0. We see from 
Eq. [5] that this represents an X error with probability 
p x = ^. From Eq. [3j 



and the total map for this noise is 

n bh [i, i,x,x\ = g = g"+g IX = g u +g XI = [1, 1, x, x 
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It is evident that this code is not very useful at correct- 
ing bit flip errors, since it leaves the bit flip channel un- 
changed. 

If the code is concatenated with itself, the resulting 
map on the bit flip channel is iW 2 o Q °* a [1, 1, x, x] — 
[1, 1, x, x], which is again ineffective at correcting bit flip 
errors. The problem is that the recovery operators are 
also concatenated using this method. Let the recovery 
operators be R = {II, XI}. Then, for one of the syn- 
dromes, we apply the recovery operator XXXI instead 
of the optimal IIIX. This syndrome has no error de- 
tected in the first block, an error detected in the second 
block, and an error detected at the top level of the code. 
The block noises are A/i = Q 11 and M% = Q XI . Applying 
the default XI recovery operator at the top level of the 
code (total recovery operator of XXXI) gives a logical 
noise of 



1, 



G = -[II - ZZ,XX - YY,XY - YX,IZ - ZI] 

l,l + x 2 l-x 2 n l + x 2 l-x 2 

= ^ — ^ ^ xO, — xQ, 

2 L 2 2 '22 

1 + x 2 1 - x 2 I + x 2 1 - x 2 , 
— ~ — — x — - — , — - — — x- 



1-x 4 



[1,1, 



2x 



2 

2.x 
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For this quasi-channel, px > Pi and py — Pz = (see 
Eq. O, and so the optimal recovery operator differs by 
encoded X from this. Since X = XX, this results in a 
recovery operator of IX instead of XI at the top level of 
the code (total recovery operator IIIX). This gives 
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The difference A = Q 11 — Q XI between these is 

1_~2 



A 



[0, 0, YX - XY, ZI - IZ] = [0, 0, x- 



x- I 
— ,x— 



Since this was the only suboptimal case, we optimize the 
total 4 qubit code by adding A to the unoptimized chan- 
nel map [I, l,x,x], yielding the optimized map 

fi 6 * [1, 1, x, x] = [I , I, Z -x - \x\ \x - l -x% 



Since x — ► |a; — ^x 3 converges to 1 for x > 0, iterating 
this map now yields a useful code for correcting X errors. 



IX. ERROR CORRECTING CODE 
THRESHOLDS 

In general, the thresholds of correctable noise for a 
code under optimized recovery operators can be esti- 
mated by analyzing for which physical noise the entropy 
of the resulting logical noise is equal to some specific value 
near the threshold. Under repeated concatenation, these 



TABLE III: Critical values for depolarizing channel 

{px,PY,Pz) = (P,P,p)- 
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d. 2979584370 


o. 2671458070 
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6.29850925% 


6.268(8)% 
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6.299(0)% 


6.269(6)% 


5 


6.299(3)% 


6.270(0)% 


6 


6.299(5)% 


6.270(3)% 


7 


6.299(6)% 


6.270(3)% 


oo 


6.299(6)% 


6.270(3)% 


Unoptimized 


4.58758548% 


3.22981197% 


TABLE IV: Critical values for a channel with independent 


probabilities p of bit flip and phase flips (p 


2 2 2\ 
~P ,P ,P~P ) 


Level 


[[5,1,3]] 


[[7,1,3]] 





lf.00278644% 


11.00278644% 


1 


10.94668310% 


10.94286393% 


2 


10.94728109% 


10.95683308% 


3 


10.949(1)% 


10.960(0)% 


4 


10.949(9)% 


10.961(5)% 


5 


10.950(4)% 


10.962(3)% 


6 


10.950(7)% 


10.962(7)% 


7 


10.950(8)% 


10.962(9)% 


oo 


10.951% 


10.963% 


Unoptimized 


7.14780025% 


6.45962393% 



critical values converge to a threshold value. Generally, 
the entropy of the physical noise at the threshold is not 
1. Noise below or above the threshold tends toward a 
logical entropy of (no noise) or 2 (totally depolarizing) . 
At the threshold, the entropy at each level of the code 
varies a lot less than the average encoded error. This al- 
lows for more efficient calculations of the thresholds than 
using the error rate We disprove the conjecture [f| 
that depolarizing noise can be corrected up to an entropy 
of 1 for the [[5,1,3]] code. The critical values are shown 
in Tables ITTT1 and HVl These results indicate that the first 
two levels, which were calculated exactly, provide a good 
estimate of the error thresholds. Higher levels were eval- 
uated with Monte Carlo calculations; for more details see 
6]. In Tab. Hill it can be seen tha,t for the depolariz- 
ing noise, the [[5, 1,3]] code has a threshold (level oo) of 
p = 6.299(6)%, which is less than the physical noise with 
an entropy of f (level 0) of p = 6.30965616%, disproving 
the conjecture of [5J that the threshold rate is equal to 
the latter. 

The value of p for which the channel Shannon entropy 
is f was calculated for various levels of concatenation 
for two different types of noise: the depolarizing chan- 
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nel (Tab. and independent bit and phase flips (Tab. 
IIVP . Tables IIIII and IIVI compare our thresholds to the 
unoptimized thresholds obtained from the channel maps 
in Refs. 0, 0|. The p values converge to the channel 
threshold for correctable noise at the oo level. The en- 
tropy at the thresholds for the [[5,1,3]] code [7| ranges 
from 0.993 (for (px,PY,Pz) — ip,0,p) noise) to over 1. 
For the [[7, 1,3]] code Q, both the phase flip noise (0, 0,p) 
and the channel (p — p 2 ,p 2 ,p — p 2 ) have a critical value of 
p = 10.963%, yielding minimum and maximum entropy 
of the thresholds as 0.4988 and 0.9976, respectively. 

The results indicate that this adaptive concatenation 
approach greatly outperforms the concatenated channel 
map method. Tab. IIIII shows that, for depolarizing 
noise, the value of p for the [[7,1,3]] code is increased 
from p = 3.2298% to p = 6.270(3)%. The [[7, 1,3]] code 



is particularly useful for quantum fault tolerance. It is 
therefore of interest to apply this adaptive concatenation 
method and similar entropy-based methods of threshold 
calculations to evaluate fault tolerance thresholds @. 

For a discussion of which noise is correctable under the 
best code for correcting that noise, see 0. 
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